load("JointQOIsHE5000.Rdata")

delta.joint  <- QOIs[,2]
rm(QOIs)

load("DirectQOIsHE5000.Rdata")

delta.direct  <- QOIs[,2]
rm(QOIs)

load("SSTQOIs5000.Rdata")

delta.sst  <- QOIs[,2]
rm(QOIs)



BP<-boxplot(cbind(delta.direct,delta.sst,delta.joint),names=c("DQ","SST","Joint"),plot=F)
bxp(BP,horizontal=T,main=expression(hat(Delta)["95%,5%"]),xlab=expression(Delta["95%,5%"]=="-0.23"),cex.main=1.4,cex.lab=1.4,cex.axis=1.4,medlwd=6,boxlwd=2,whisklwd=2,staplelwd=2)
lines(c(-0.225511,-0.225511),c(0.43,1.54),col="gray80",lwd=7)
lines(c(-0.225511,-0.225511),c(2.455,2.54),col="gray80",lwd=7)
lines(c(-0.225511,-0.225511),c(3.46,3.56),col="gray80",lwd=7)



